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Abstract 

A simple model is presented to calculate the potential of mean force between a polyion and a 
multivalent counterion inside a polyelectrolyte solution. We find that under certain conditions the 
electrostatic interactions can lead to a strong attraction between the polyions and the multivalent 
counterions, favoring formation of overcharged polyion-counterion complexes. It is found that small 
concentrations of salt enhance the overcharging while an excessive amount of salt hinders the charge 
reversal. The kinetic limitations to overcharging are also examined. 



I. INTRODUCTION 



It is our pleasure to contribute this paper to the special issue of Molecular Physics dedi- 
cated to celebrate Ben Widom's outstanding contributions to Physical Chemistry and Sta- 
tistical Mechanics. Ben's work is characterized by a profound physical insight, combined 
with an ability to abstract the most complex physical phenomena into a simple model. 
From scaling and criticality^ to microemulsions^ and the hydrophobic effect^, Ben's sagacity 
has opened new frontiers of Physical Chemistry. While it is impossible to compete with 
Ben's intuition, one can at least try to follow his example. In this paper we will, therefore, 
study a simple model of interaction between a polyion and multivalent counterions inside a 
polyelectrolyte solution. 

Thermodynamic systems in which long range Coulomb interactions play the dominant 
role pose an outstanding challenge to Physical ChemistryA. Even such basic question as the 
possible existence of a liquid-gas phase separation in a restricted primitive model has been 
positively settled only quite recently^. Even so, the order of this transition still remains 
a source of an outstanding debate^. For strongly asymmetric electrolytes such as aqueous 
colloidal suspensions, even the existence of a liquid-liquid phase separation continues to be 
controversial^i^iSiiSiii. 

When aqueous colloidal suspensions or polyelectrolyte solutions contain multivalent coun- 
terions other curious phenomena appear. For example, it is found that for sufficiently small 
separations two like-charged polyions can attract one anothe r'^^1^^1^^1^^1^^1^'^1^^1^^ . If an ex- 
ternal electric field is applied to such a suspension the electrophoretic mobility of colloidal 
particles is often found to be reversed, so that the particles move in the direction opposite to 
the one expected based purely on their chemical charge^i^SiSLiSiS^. Both of these phenomena 
are a consequence of strong electrostatic coupling between the polyions and the counterions. 

The counterions inside the suspension can be divided into two categories: those which 
are associated (condensed) to the colloidal particle and those which are free. The condensed 
counterions contribute to the effective, renormalized, charge of the polyion-counterion com- 
plex, while the free counterions and colons result in screening of the electrostatic interactions 
inside the suspension^. In this paper we will explore the potential of mean force between 
a rodlike polyion with n associated counterions and a counterion located at a transverse 
distance d from the polyion center. Fig. ^ 
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FIG. 1: Rodlike polyion with Z = 6 monomers separated by a distance b and a counterion located 
at X = d. 

II. THE MODEL 

Consider a rodlike polyion of Z (even) monomers, each carrying a charge —q, inside an 
aqueous suspension containing multivalent counterions and salt. The monomers are located 
uniformly with separation b along the rod. Strong electrostatic coupling between the polyion 
and the counterions results in a condensation of n a-valent counterions onto the polyion. The 
condensed counterions are free to hop between the monomers of the polyioi>i^. If a monomer 
has an associated counterion, its charge is renormalized to {a — l)q. The free, uncondensed, 
counterions and colons screen the electrostatic interactions, changing the potential between 
the two charges qi and q2 from the Coulomb to the Debye-Hiickel-^ form 



where e is the dielectric constant of the solvent and k is the inverse Debye length. The 
question that we would like to address in this paper is what is the potential of mean force 
between the polyion-counterion complex containing n condensed a-ions and an additional 
a-valent counterion located transversely at distance d from the polyion center, see Fig. 1. 

To proceed, we assign to each monomer i a lattice-gas variable Uj, such that cij is equal 
to 1 if a counterion is condensed onto cite i and otherwise. For a given configuration {a}, 
the interaction Hamiltonian between the complex and a counterion located at a transverse 
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distance d from its center is 
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where 
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It is convenient to define the reduced distance between the polyion and the counterion 
X = d/b, the reduced inverse Debye length k = nb, and the Manning paramete r^^i^^ as ^ = 
q'^/eksTb. In terms of these adimensional variables the reduced Hamiltonian, H = /37i/^, 
becomes 



H = - 1) 
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{2i-l- Zy + 4x2 
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The partition function is a trace over all possible distributions of n condensed counterions 
among the Z polyion sites. There is a total of 

Z\ 



Nr. 



(Z — n)\n\ 

such configurations. The partition function is then 

where the sum is over the Nc configurations {a} which obey the constraint J2i=i = 
denoted by the prime. It is convenient to order the terms in the Hamiltonian by the distances 
between the pair of interacting charges. This results in 
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If we now define the Boltzmann factors 



Xj 



exp 



-^exp(-A:./(2j - 1 - Zy + Ax^) 



and 



Vj = exp 



(2j-l-Z)2 + 4x2 
^exp(-A;j) 



J 



the contribution of each configuration to the partition function will be a product of these 
factors raised to exponents which are polynomials in a, that is 



z~i 



Nc Z/2 
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The polynomials Vij = —aija + bija'^ and Uij = Cij — dija + Cija^, have integer non-negative 
coefficients. The advantage of the simple model constructed above is that for not too large 
values of Z and n the partition function can be evaluated exactly with a help of a computer. 

The potential of mean force (measured in units of /eh) between a polyion-counterion 
complex and an a-ion located at x is, 



0(4, = --lu' 



(5) 



i g(oo) 

The potential is normalized so that ^(cxd) = 0. 

The computer code which generates the partition function for given values of Z and n 
determines the set of integer coefficients of the polynomials defined following the Eq.(@]). 
Each set of polynomial coefficients may correspond to more than one internal configuration 
of the polyion, so that the degeneracy must also be taken into account. All the data is stored 
on the computer and used to perform a floating point calculation of the free energy. 



III. RESULTS AND DISCUSSION 

In Fig. 121 the potential of mean force between various complexes and an a-ion is plotted. 
The complexes are composed of a polyion of charge — lOg and n associated divalent counte- 
rions. Notice that for n = 5 (neutral complex) the potential is a monotonically increasing 
function of x, so that the sixth counterion is always attracted to the complex. 
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FIG. 2: Potential of mean force as a function of x for Z = 10, ^ = 1, /c = 0, and a = 2. In upward 
order, the curves correspond to n = 5, 6, 7, 8, 9 condensed counterions. 

For an overcharged complex with n = 6 condensed counterions, the potential of mean force 
develops a barrier. At large distances the seventh counterion is repelled from the complex, 
while at short distances it is attracted to it. The minimum of the free energy, however, 
is reached when the seventh counterion is located at x = 0. The potential of mean force, 
therefore, favors counterion condensation. The size of the barrier increases with n and the 
minimum at x = becomes metastable for n = 8. For n = 9 the potential is a monotonically 
decreasing function of x, and the tenth counterion is always repelled from the complex. We 
next study the dependence of the depth of the potential well and the height of the barrier 
on the parameters of the model. 

In Fig. we saw that when the complex is overcharged n > Z/a, the potential can have 
two minima, one located at x = oo and another at x = 0. Which one of the two minima is 
the global one is determined by the sign of 0(0). Figs. El and 0] show the behavior of 0(0) 
as a function of ^ and k. When 0(0) < the position at x = is the absolute minimum, 
while when 0(0) > 0, x = is at most metastable. We should note, however, that the 
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FIG. 3: The potential of mean force at x = as a function of ^ for some values of k. The curves 
are for Z = 10, n = 4, and a = 3. 

present discussion is not sufficient to define the absolute number of condensed counterions. 
For a counterion to be condensed the depth of the potential well must be sufficiently large, 
compared to the thermal energy fcsT, to prevent its rapid escape from the polyion surface. 
At the level of the present discussion this criterion is arbitrary. Thus, in this paper we will 
not consider the absolute number of condensed counterions but only the conditions which 
favor or disfavor the counterions condensation. From Figs. ElandlU we see that for a polyion 
of Z = 10 and n = 4 condensed trivalent counterions the minimum at x = is the global 
one for the parameters plotted. Approach of an additional fifth counterion to this already 
overcharged complex is, therefore, energetically favorable. 

The depth of the global minimum |0(O)| is a monotonically increasing function of the 
Manning parameter, see Fig. 01 The dependence on the salt concentration, however, is 
not monotonic. From Fig. |3] we see that small concentrations of salt favor counterion 
condensation, i.e. 0(0) becomes more negative for small k. Larger concentrations of salt, 
however, have a destabilizing effect on the counterion condensation. 
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FIG. 4: 0(0) as a function of k for some values of The curves are for Z = 10, n = 4, and a = 3. 

This is even clearer for complexes composed of a polyion with Z = 10 and n = 5 condensed 
trivalent counterions. Fig. |S1 shows that the position of the free energy minimum is a 
non-trivial function of salt concentration. Depending on the Manning parameter ^ and the 
concentration of salt k, association of an additional, sixth, counterion can be either favored 
or disfavored. On the other hand, for Z = 10, n = 6, and a = 3, 0(0) is always positive so 
that a complex with n = 7 condensed counterions can be at most metastable. 

We next explore the dependence of the barrier height 0(xm), where Xm is the position of 
the maximum of the potential of mean force see Fig. |2l on the parameters of the model. 
In Fig. ini (pixm) is depicted as a function of the Manning parameter ^ for polyion of size 
Z = 10 with n = 4 associated counterions. We see that the barrier height diminishes with 
the increase of ^ and the amount of salt inside the suspension. 

To explore the dependence of the barrier height on the size of the polyion Z, in Fig. [7|we plot 
(p{xm) as a function of Z for complexes composed of a polyion and n* condensed trivalent 
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FIG. 5: 0(0) as a function of k for some values of ^. The curves are for Z = 10, n = 5, and a = 3. 
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FIG. 6: The barrier height 4>{xm) as a function of ^ for some values of k. The curves are for Z = 10, 
n = 4, and a = 3. 
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FIG. 7: The barrier for polyions with ^ = 1 and n* condensed trivalent counterions, such that 
<Pn* (0) = 0, as a function of Z. 

counterions, such that 0n*(O) = 0. While in the absence of salt the barrier height shows 
a significant dependence on the polyion size, at finite salt concentration this dependence 
weakens and (f){xm) seems to saturate when the polyion size is significantly larger than the 
Debye length. For large Z and small concentration of electrolyte, however, the kinetic barrier 
can be many ksT, providing a significant limitation to overcharging^^. 



Charge reversal is a consequence of strong positional correlations between the counterions. 
These correlations are induced by the electrostatic repulsion between the particles. Thus, 
we expect that both the barrier height and the relative depth of the absolute minimum will 
be strongly dependent on the counterion valence. In Figs. |8land|9lwe show the dependence 
of the barrier height and the depth of the potential well on the valence of the counterions. 
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FIG. 8: The barrier for complexes composed of Z = 20 and n condensed counterions, such that 
na = 24, as a function of the counterion valence a. 

Although all the overcharged complexes depicted in Figs. |H1 and IHl have the same net charge 
4g, the depth of the potential well and the height of the kinetic barrier depend on a. As 
expected, larger counterion charge leads to stronger positional correlations and favors the 
counterion condensation and the charge reversal (0(0) becomes more negative with increasing 
a). The barrier height, however, once again shows a nontrivial dependence on the salt 
concentration. For small amounts of salt and large Z, increased counterion valence leads to 
larger kinetic barriers. 

IV. CONCLUSIONS 

We have studied the potential of mean force between a polyion and an a-valent counterion 
inside a polyelectrolyte solution containing multivalent counterions and monovalent salt. 
The model is sufficiently simple that the partition function can be calculated exactly. It is 
found that for an overcharged polyion the potential of mean force can have two minima, one 
located at x = and another x = oo. Which one of the minima is the global one depends on 
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FIG. 9: (p{0) for complexes composed of Z = 20 and n condensed counterions, such that na = 24, 
as a function of the counterion valence a. 

the charge density of the polyion and the amount of salt inside the suspension. When the 
global minimum is at x = 0, a counterion from the bulk finds it energetically favorable to 
approach the polyion surface. To reach x = 0, however, the counterion must overcome a free 
energy barrier. For small salt concentrations, this barrier can be sufficiently large to provide 
a kinetic limitation to the extent of charge reversal. Furthermore, even if the counterion 
reaches a; = 0, whether or not it will become condensed will depend on the depth of the 
potential well. Counterion condensation will occur only if (j){xm) — 0(0) ^1/^. Otherwise, 
the thermal fluctuations will lead to a fast escape of the counterion from the x = minimum. 

For suspensions containing rodlike polyelectrolytes and the multivalent counterions micro 
phase separation is observed under certain condition o^^'^^ . The polyions aggregate forming 
bundles with a well defined crossectional area. It has been argued that bundle formation is 
an activated process and the size of the bundles is kinetically controlle d^°i^^i^^ . It should then 
be quite interesting to explore the dependence of the barrier height on the concentration of 
monovalent electrolyte using a theory similar to the one presented above. 

To conclude, the extent of the charge reversal is strongly dependent on the amount of 



12 



monovalent salt present in the suspension. Small concentrations of salt will enhance the 
overcharging while an excessive amount of salt will hinder the charge reversal. Furthermore, 
even if the minimum of the free energy corresponds to an overcharged state, we find that 
depending on the polyion charge density and the amount of salt in the suspension, there can 
be significant kinetic limitations to the overcharging. 

This work was supported in part by the Brazilian agencies CNPq and PAPER J. JFS 
acknowledges funding by project Pronex-CNPq-FAPERJ/171. 168-2003. 
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